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ABSTRACT 


In its motion response to a disturbance, a ship can be closely ap- 
proximated as a second order system. The solution of these second order 
equations is used to determine the motion stability of the ship. By ap-= 
plying a step disturbance to a ship model and measuring the response, 
the linearity can be used to find the response to any disturbance, 


In this work a step disturbance was applied to a Mariner class 
hullform and the response measured. From this response the response to 
regular sinusoidal excitations was computed. The response to the sinu= 
soidal inputs was then used to determine the coefficients of virtual 
inertia and damping for the Mariner in pitch and heave, 


The derived results at zero speed were compared with values com= 
puted by strip theory. 
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CHAPTER I 


INTRODUCTION 


The equations of ship motion developed by Abkowitz (1) are a 
linearized set of second order partial differential equations. In the 


vertical plane the equations are 


(m-Z..)z - Z.2-Z 2 - (Z,+nX,)6- - (Z+uZ.,)8 - (2,+u,2, )6 =a ACr a) (1) 


(I, -M)8 - (Mstu M, )@ - (MtuM, ) 6 - (M,tmiX  )z - M,Z -M2 = M(ex) (2) 


z = displacement in the vertical plane 
zZ = velocity in the vertical plane 

Z © acceleration in the vertical plane 
6 = angular displacement in pitch 

6 = angular velocity in pitch 

6 = angular acceleration in pitch 
Uo = forward velocity of the ship 

m = mass of the ship 
Z-- = added mass coefficient 
Z- © damping coefficient 
Z = static restoring force coefficient 
X. = distance from origin to the center of gravity(equal to zero) 


Z- ® coupled virtual mass 


S 
Ze = rotary force derivative coefficient 
Ze = static force derivative coefficient 
De = moment of inertia 
Ma = added moment of inertia coefficient 
Me = damping coefficient 


=i 





Mg = static restoring moment coefficient 
M;- = coupled virtual moment of inertia coefficient 
M. = rotary moment derivative coefficient 
M © static moment derivative coefficient 


These equations show a coupling effect between pitch and heave, 
However, by decoupling these equations and neglecting second order 


effects, these equations can be reduced to 


4 


Gee Ze. )2 - Ze 2 -Z25 Z(ex) (3) 


-_ 


ce : Mix)8 M0 - M,6 = M(ex) (4) 


The equations are now in a fom that is readily solvable. The 
solution of these uncoupled equations is highly significant in deter- 
mining the motion characteristics of a ship. These characteristics are 
important in determining stability of motion as well as factors affecte 
ing the comfort and safety of personnel embarked on the ship, Before 
these equations can be solved, however, it is first necessary to deter~ 
mine the hydrodynamic coefficients of added inertia and damping, 

The added inertia is a term that results from the fact that a ship 
moves in a dense medium and therefore must displace a heavy mass of 
liquid whenever it moves in any plane, This movement of the mass of 
liquid has the effect of essentially increasing the mass of the ship, 
The damping is a term associated with energy dissipation and not vis~ 


cous effects, As a surface ship operates at the interface between water 


wo? ws 








and air, any vertical motion results in the creation of waves. The 
energy lost in wave creation acts as a damping force to reduce the os- 
cillations of the ship, . 

Many methods, both theoretical and experimental, have been de- 
veloped to determine these coefficients, One such method was developed 
by Narita (11) at the haseachivsetts Institute of Technology in 1962, 

It consisted of applying a known moment step disturbance to a hull form 
and measuring the transient response in pitch, while restraining the 
model from motion in heave, The resultant motion of the model was a 
transient wave of a aenapine Sinusoid form leading to the step. When 
this transient response was obtained, it was used to compute the theo= 
retical response of the model to any known input, This makes use of a 
property of a linear system, that when its response to a step or im= 
pulse input is known, its response to any other input may be found by 
use of the convolution integral, As a ship is closely approximated by a 
linear second order system, this principle is applicable. 

This property can be used to apply known regular sinusoidal inputs 
and determine the ship's response. This in essence produces the same 
results as could be obtained by making forced oscillation tests, The 
output of the regular sinusoidal inputs can then be used to calculate 
added mass and damping. 

Narita (11) computed added inertia and damping for the pitch 
motion only. In this work the procedure was extended to measure the co= 
efficients for pitch and heave, Two series of tests were run, one in 


which the model could pitch with heave restrained and one in which the 





model could heave with pitch restrained, 

The authors know of no previous tests conducted on the Mariner 
hull form for the purpose of determining coefficients of damping and 
added mass. However, a theoretical prediction for the zero-speed con- 
dition may be approximated. Porter (1) has reported that the conformal 


transformation of the unit circle represented by the equation 
N 


ee li (5) 


n=0 
where z is the plane of the actual ship section and ¢ is the plane of 
the unit circle, can en a given station of the ship. N is se» 
lected and the coefficients ao») are found which give the best fit to 
the actual section, Lewis (10) used this transform with N = 1 and 
Landweber and Macagno (9) extended it to N= 2, In the present work, 
a method was used which can give a least mean square error approximation 
to the actual ship section if the number of coefficients is selected. 
The least number of coefficients which gave a reasonable fit for each 
section was used, 

After the coefficients of the transform were determined, Porter's 
method was used to obtain coefficients of damping and added mass for 
each station at several frequencies and the results were numerically ine 
tegrated throughout the length of the model by using strip theory to 


obtain added mass and damping for the model, 





CHAPTER ITI 
PROCEDURE 


2.1 Model Characteristics 
The tests were conducted using a model of the Mariner class hull 


form, which has the characteristics shown in Table I, 


Table I 
CHARACTERISTICS OF THE MODEL 


od 


Length between perpendiculars _ 5.28 feet 

Beam 76 feet 

Draft o2( feet 

Displacement in fresh water 2. pounds 

Block coefficient 06125 

Radius of gyration 1.26 feet 

Static restoring moment (pitch) 158.8) foot - pounds/radian 
Static restoring force (heave) 148,80 pounds/foot 

Material of model Fiberglass, reinforced plastic 


2e2 Test Conditions 

The tests were conducted at the Massachusetts Institute of Teche 
nology ship model towing tank. The tests were divided into two series, 
pitch with no heave and heave with no pitch, Each series consisted of 
tests conducted at Froude numbers of 0.0, 0.1, 0.2, 0.3. Further, the 
tests conducted at zero Froude number were conducted parallel to and 


perpendicular to the centerline of the towing tank in order to investigate 


52 


the effects of wall reflections, In the pitch -no heave, configuration, 
the ballast weights were shifted in order to determine the effect of 


changing the radius of gyration on the computed coefficients. 


2.3 Step Disturbance Application 
2e3el General 

As the tests were to be conducted in a liquid medium, it was de~ 
cided that the high potential devices, as recommended by Narita (11), 
posed a danger to personnel working with the equipment and to other 
electronic equipment in the event of a wiring mistake or a faulty lead. 
For this reason low current and voltage devices were investigated, and 
it was decided that the best device might be a small electromagnet. The 
magnet chosen for these tests was a simple iron core wound magnet. For 
the purposes of these tests, the magnet required less than two amperes 
of current and less than two volts of potential. These values were not 
only safe to personnel in the event of a casualty, but were easily ob= 
tainable with standard direct current power supplies, Initially bate 
teries were employed, but as the tests took a long time, because of de- 
lays in waiting for the tank to settle to a dead calm, the batteries 
rapidly lost their potential. In place of the batteries a Sorenson 


Model QR 18-1.5A direct current power supply was used, 


2.3.2 Pitch-—No Heave Configuration 
To apply the pitch disturbance, the model was mounted as shown in 
Figure I. A small metal plate was mounted on the bow and on the ster 


of the model, <A millimeter bar was used to mount the magnet. This 
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Figure I 


PITCH-NO HEAVE CONFIGURATION 
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enabled the operators to have greater ease in setting the disturbance at 
very close to one and a half degrees, As mentioned above, ‘the current 
was supplied to the magnet from a Sorenson Model QR 18-1,5A direct cur- 
rent power supply. The electric wiring arrangements were as shown in 
Figure II, By releasing the current to the magnet, the ship will let go 
under a step-like displacement, following a decaying sinusoid coming to 
rest at a zero angle of pitch, 

As the exact moment is desired at which the ship actually reacts 
to the disturbance, a second circuit was constructed as shown in 
Figure III, When the ship hee to react, the magnet parts from the 
steel plate and produces a step reproduction on the recorder, The be- 
ginning point of this step can be used to determine the starting point 
of the model's response and thus increase the accuracy of calculation 


based on the response of the model, 


2.3.3 Heave—No Pitch Configuration 

In constructing the heave arrangements, it was necessary to devise 
a method of lifting the entire ship, It was also desirable to not ine 
troduce any forces tending to cause pitch, which would thus put undue 
strain on the heave rod, For this reason it was decided to mount a 
small metal plug directly on the heave rod and to lift the model by 
connecting the magnet to the metal plug, As the heave rod is mounted 
approximately at the longitudinal center of floatation of the model, es- 
sentially no pitch moment is introduced. The magnet was fastened to the 
carriage by two aluminum rods as shown in Figure IV. The circuitry is 


exactly the same as for the pitch configuration. 
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Figure II 
ELECTRIC WIRING ARRANGEMENT 
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Figure IV 
HEAVE-NO PITCH CONFIGURATION 
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In both heave and pitch testing the magnet and metal connecting 
plate were mounted rigidly, so there was no concern for catching a falle 
ing weight. Then to reconnect, it was only necessary to close the 
switch and lift the model so that the plate and the magnet came in cone 
tact, at which time the magnet would support the weight of the model, 
Then it was ah iiclvouayate wait until the tank settled to an absolute 


calm before beginning the next run. 
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2.45 Measurement of the Step Response 
2euet Piteh 

In order to increase the validity of linearity of motion theory, 
an initial displacement of only one and a half degrees was used. To get 
the maximum accuracy in measuring this small variation, two methods were 
- investigated, use of a linearsyn and use of a pitch bearing, To use a 
linearsyn, maximum excursion of motion is desired, This would be at the 
ends of the model. However, the motion here is circular enough to cause 
non=linear motion of the linearsyn and create inaccuracies, Also, the 
vibrations near the ends of the model were found to be the greatest, and 
this too would have an adverse effect on measurements made by a linear~ 
syne 

A pitch bearing located right at the longitudinal center of 
floatation has the advantage that it can measure angular motion right 
at the point about which pitch takes place, The pitch bearing is a 
rotary motion displacement pickoff of a differential transformer type 
and is highly sensitive to any rotary displacement, Therefore, it was 


decided to employ the pitch bearing for the pitch measurements, 


ye 





2.4.2 Heave 

To measure heave, the linearsyn was perfectly suited, It is a 
linear motion displacement pickoff of the differential transformer type, 
It consists of a coil assembly and movable magnetic core and is capable 
of detecting motions as small as 0.000001 of an inch between the core and 
coil assembly, Mounted directly onto the heave rod assembly, the 
linearsyn gives an accurate reading as to the movement of the model in 
the heave plane. j 

To increase the validity of model Linearity in heave motion, a 
disturbance of only four-tenths of an inch was applied to the model, 
This proved sufficient to obtain several well-defined cycles of oscil- 


lation, 


2.5 Test Procedure 

After the model was mounted on the carriage and circuitry arranged, 
the magnet was connected to either displace the model in pitch or heave, 
depending on which tests were being conducted, Tests were then carried 
out in two groups, pitch and heave, at Froude Numbers of 0.0, O.1, 0.2, 
0.3. 

With the model displaced it was necessary to wait until the tank 
had settled to an absolute calm before beginning a test. Then the model 
was brought up to speed and the switch opened, stopping power to the 
magnet, At that instant the ship would fall toward zero in a decaying 
sinusoidal motion, characteristic of a second order system. This motion, 
as well as the exact instant at which the magnet parted from the metal 


plate, was recorded on both the Sanborn recording paper and on magnetic 


joe 





tape for later analysis, To begin the next run, the model had only to 
be towed back to the opposite end of the tank, the magnet re~engaged 
with the metal plate, and then allowed to remain suspended until the 
tank settled, 

For the zero speed runs the model was placed at the center of the 
tank, first parallel to and then perpendicular to the centerline, Here 
the procedure was the same as before except that there was no forward 
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motion of the model, 


a4 3. 


CHAPTER III 


ANALYSIS 


3.1 Background 
The analysis procedure used in this work was essentially that used 
by Kerwin and Narita (5). A step input is applied to a system, and re- 


sponse is represented by 
P(t) = = (1 + R(t)) (6) 


where C is either the stateeP eau risic force in heave or the static ree 
storing moment in pitch, R(t) represents the transient portion of the 
response and 1/C is the final value the response seeks after the tran- 
sient has died out, 

To get the response to a regular sinusoidal input, the convolution 


integral is used. 


& 


v(t) = fv (r) H(b = 1) ar (1) 


H(t) represents the impulse response of the system, However, as the step 


response, S(t), is known, it may be used with the input in the form 





© daVv.(r) 
V(t) 38 i S(t - +) dt (8) 


any 
and the same results will be obtained. From this, the response of the 
model to sinusoidal inputs in either heave or pitch may be obtained, 


To obtain the damping and virtual mass coefficients, we represent 


ea Ne 





our second order equations as follows: 


2 
qd 6 ade : 
ie pap | “pes M. sinwt (9) 
at 
2 
a 2 az . 
“HE até eS: Br (10) 


The solution of these equations will be of the form 


a sinwt + b coswt (11) 


from which 





2 2 
—_ Cy - M9,/ (a, +b, ) is 
P 2. 
W 
-M b 
_ one 
i (aie Spe) (13) 
ee P 
2 2 
Clea Leama be | 
H O H lal 
ty Ey Be” (a) 
W 
-Z2 0D 
B= wollte seers (15) 
ua ae) 
H ig 


302 Static Coefficients 

Equation (6) indicates that the step response tends to 1/C as the 
transient dies out. In the initial simplification of the equations, it 
was assumed that this was a static term and not speed dependent. Ger= 
ritsma (2) shows that this term is speed dependent and that this depend-= 


ency is further a function of the hull form, While the variation in 


aie. 





speed is not great at very low Froude numbers, it does approach approxi- 
mately 10 percent at a Froude ee of three-tenths, This is still not a 
large effect, but it may exert a measurable influence on the values of 
added mass and damping calculated ae the various Froude numbers, 

Accordingly, tests were conducted to determine the change in the 
static restoring coefficients with speed, In pitch a known moment was 
applied to the ship and the pitch angle was measured, As the pitch 
velocity and acceleration were *both zero, the static restoring coeffi- 
cient could be found by dividing the moment by the resultant angle. The 
static restoring moment at zero speed was computed to be 158.8 ft.-lb./ 
radian, This same procedure was repeated at Froude numbers of O.1, 0.2, 
and 0.3, and a variation of less than 1 percent was observed . each 
speed. Therefore, the static coefficient at zero speed was used for the 
analysis at each of the Froude numbers, 

A known heaving force was then applied to the model, and heave was 
measured at Froude numbers of 0.0, 0.1, 0.2, 0.3. In this case the re- 


sults followed the form of those found by Gerritsma as shown in Table II, 


Table il 


STATIC RESTORING CORFFICIENTS IN HEAVE 


ee ees 2 
0 148,8 lbs/ft 
i 11.0 lbs/ft 
ae 1h2.8 1lbs/ft 
oo 139.2 lbs/ft 
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The variation is not large, but is significant enough to warrant 


uSing the value found at its corresponding Froude number, 


323 Digitizing Rate 
The data as obtained from the experiments was recorded in analog 


form on magnetic tape. As the convolution integral was performed by 
series methods, it was necessary to have the data in digital form, Pre- 
vious experiments conducted by Kerwin and Narita (5) showed the results 
of added mass and damping to be highly dependent on the digitizing rate 
used, with damping being the more critical. From their tests of a 

known second order system, it was found that a digitizing rate of one 
kilocycle per second gave exact results for added mass and to within 
three--tenths of 1 per cent for damping, 

Accordingly, a digitizing rate of one-thousandth of a second was 
used, To digitize the data at this rate, it was necessary to record a 
one kilocycle per second square wave on the tape, as the process used 
with the IBM 709) computer reads the value of data everytime there is a 
positive to negative vertical change in the digitizing signal, 

As mentioned earlier, the response of the model was recorded as 
well as a signal that showed the instant the magnet parted. These sig- 
nals were recorded first, and the square wave was recorded afterward, 
The square wave was recorded insuring that it started prior to the in=- 
stant of the magnet parting signal. By this procedure and obtaining the 
digitized data for both the response data channel and the parting signal, 
the parting signal data could be scanned, and the exact instant of mage 


net parting could be found, On the other channel, this was the starting 
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point for which the response signal data was used in the analysis pro- 
gram. 

As mentioned above, to get maximum accuracy, a digitizing rate of 
a thousand points per second was used to obtain the best representation 
of the model response data. In convolving this data with the input 
sinusoid, a series method was used, This involved shifting the sinu- 
soid in reference to the response data. The response to the sinusoid 
was computed as a series, and the final value summed. In shifting the 
waveform, it was found that the fine spacing of .00l1 seconds was not 
required. The analysis program was altered to determine what interval 
could be used and still produce the same results that one-thousandth of 
a second produced, The main purpose of this alteration was to speed 
up the analysis program and reduce computer time used in obtaining re- 
sults. Table III shows the results of increasing the convolution 
interval at several frequencies or corresponding periods of excitation, 

These values were obtained using ideal data for the second order 
— 26 “sg e575 = + 1776 (16) 

From Table III it can be seen that increasing the convolution ine 
terval actually improved the accuracy of results at all except the high- 
est frequency, period equal to 0,6) seconds, In this case the added in-= 
ertia stayed approximately the same, but the damping became worse by 
seven=thousandths, The over-all results did show that larger convolution 


intervals gave accurate results. Accordingly, larger intervals were used 
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in the analysis, depending on the input excitation frequency. 


Table aT, 


EFFECT OF CONVOLUTION INTERVAL ON COEFFICIENTS 


Period Convolution Interval Added Inertia Damping 
(Seconds) ~ (Seconds) (Slug Feet2) (Slug Feet?) 
1.28 sO | 2507 “a ; 581 
1.28 e010 2.505 72550 
1.28 O10 2503 72549 
1.28 2020 2503 7.549 
a2 s0B2" 2.502 75446 
1.28 «00 2.502 755 
2.56 s0OL Zo sib 7658 
2.56 e010 2.506 7-602 
2 56 .016 2.506 7.601 
2.56 .020 2.506 7.600 
6h 001 2.198 70529 
6h, 010 2499 70536 
6h 2020 2499 70536 


3.4 Theoretical Analysis 

In the present problem offsets for 21 stations of the Mariner hull 
were read off a +-inch-scale body plan, and by iteration on the IBM 709), 
computer, the coefficients ao,4) which gave the best fit to the actual 
sections were determined, No more than five such coefficients for any 
one station were used because the physical significance of coefficients 


beyond Ag becomes unclear, The values of the coefficients giving the 


=19% 





best fit for each section are tabulated in Table IV, 

After the coefficients a5,, were determined, they were used as 
inputs to a computer program which utilizes Porter's methcd as explained 
in reference 1); and calculates coefficients of damping and added mass 
for each station. The results for each station were integrated over the 
length of the model by using strip theory to give the coefficients for 


the model, 
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Table IV 


SECTION TRANSFORMATION COEFFICIENTS 


Section Area 


os 


O 0.20 210 7407 
i oo” 2720 0.999 
2 Sus ifs 7M 0. 0.818 
3 16.15 > 27.0 0.779 
\ 22.70 27.0 0.791 
5 28.50 27.0 0,830 
6 S525 27 20 0.86) 
7 36. 30 27.0 0.91) 
8 Syewie 27.0 0.958 
9 38,00 27.0 0.982 
10 38.00 27.0 0.988 
1 38,00 27.0 0.979 
12 38,00 276 0.965 
13 38.00 A 0.932 
1h Stat 27.0 0.868 
15 36.45 ol 10 0.787 
16 3340 2160 0.700 
ii 28,25 270 0.609 
18 21.05 27.0 0.519 
19 12.20 27.0 0.390 
20 1.80 0.85 0,020 


Table IV (continued) 


00 


Station ay a3 ac a7 ag 
O -.85755187 ~.11659060 ~,0189169) +.0065186719 -~.,0003725316) 
1 ~-.660394  =-.102210 -=,006052 + ,0211))2 
2 -.h7394K1 -=,08562. +,012796 4.01726 +,006),57 
3. =.273985 = #02895 = # ,02)931 #015372 
4 = =6112573 0 -,017126 = + 019732 + ,0071,28 +,007161 
5S  #,001188  -.032628  +,019083 4.000588 +,005890 
6 +.076083 ~.05392 + 016873 + ,000628 +,00529) 
7 +.112537 -.078719 +.010882 ~ 001901 + 001656 
8 +,143558  -.103009 +,00h78) 

9 +,150)10 ~ 118791 ~,001681 ~ 002359 
10 +,150708  -,125100 -.003798 ~,002210 + 000777 
11 = #,151138 = =6-,123101 = = 00327 = 003000 
12 +,150843  -.110398  -,001079 ~ 00638 
Ls + ,1),8500 ~ 085690 + ,005),38 - 004678 
1h +*,.1h0856 -+,08439 +,015873 - .00),178 + ,000558 
15 +#,12407)  =,002855 +*,023)61 ~.001276 4.000786 
16 +,081437 +.049610 +,029596 ~ 001735 

17 -.022515 +,101750 + 01,2262 -.00002 +.005171 
18 =,220882  +,138063 +,078202 + ,001977 +,001509 
19 ~ 55295 +,110698 +,101173 + 03500) + ,020560 
20 $3192) / ~.08)215 + 05,789 =.052018 4.019595 





CHAPTER IV 


RESULTS 


4.1 Heave-No Pitch Condition 

To obtain values of added mass and damping in heave, four tests 
were made at each of the Froude numbers desired, A cee response 
curve is shown in Appendix B,. These curves are presented to show 
wall reflection effects, but they ‘also show the typical response at zero 
speed. At forward speeds greater than zero the model response curves 
were much the same, but considerable vibration was also introduced be~ 


cause of the towing carriage oscillations. 


4.1.1 Virtual Mass 

The tests conducted at each aieed were averaged to obtain the vire 
tual mass, This average value was then nondimensionalized by dividing 
it by the mass of the model. Each of these values is shown in Table V 
as a function of excitation frequency. Furthermore, each of the ob-= 
Served values was compared with the average to obtain the worst excur= 
Sion, This worst excursion from the average is expressed as a percent 
of the average in Table V. 

In most cases all points were very close together, and percentage 
errors were small, However, occasionally a derived value was at vari- 
ance with the rest of the results and its error appeared quite large. 
This took place at those values farthest from the model's natural fre= 


quency of approximately 7.6 radians per second, 


a3 








The symbols used in Table V are as follows; 
w radian frequency of excitation 
A virtual mass 
At nondimensional added mass 
g gravity 
A displacement of the model 
The average values of the virtual mass are plotted in Figure V 
for each of the four Froude riot Figure VI is a comparison of the 
zero speed values obtained with the model parallel to and perpendicular 


to the centerline of the tank, - 


4.1.2 Damping 

The same procedure was used with damping as was used with added 
mass, The results of the average values, as well as the maximum percent 
error, are shown in Table VI.. The damping values were nondimensional- 
ized by multiplying the damping Sedthea by the square root of the 
length of the model times gravity divided by the displacement of the 
model. 

Errors in damping were very small, but were generally slightly 
higher than for added mass. However, near the natural frequency of the 
model, the errors were very small, increasing as the frequency was in- 
creased or decreased, 

Additional symbols introduced in Table VI are: 

B damping 
B' nondimensionalized damping 


L length of model 
J2 hve 





The average values of damping are plotted in Figure VII for each of 
the four Froude numbers. Figure VIII is a comparison of the zero speed 
values obtained with the model parallel to and perpendicular to the cen~ 


terline of the tank, 


4.2 Pitch-No Heave Condition 

To obtain the values of the added inertia and damping coefficients 
' in pitch, six tests were conducted at each of the Froude numbers, Three 
tests were conducted with the how initially up, and three tests were con= 
ducted with the stern initially up. The response curves were the same as 
for the case of heave, in that each was a decaying sinusoid, At speeds 


other than zero there was a large amount of vibration in the response, 


paralleling the situation found in heave, 


4.2.1 Virtual Inertia 

The same analysis program was used to compute the value of virtual 
inertia in pitch as was used to compute the virtual mass in heave, The 
only changes were to change several constants and the nondimensionalizing 
factor, Virtual inertia was nondimensionalized by multiplying the value 
found by the value of gravity divided by the displacement times the 
length squared, The values of virtual inertia found at each of the four 
Froude numbers are shown in Table VII, In this table A designates the 
value of virtual inertia, Figure IX is a plot of the virtual inertia 
versus excitation frequency for each of the four Froude numbers, Fig-= 
ure X is a comparison of the zero speed condition with the model paral- 


lel to and perpendicular to the centerline of the tank, This figure 


a5 





also shows the effect on the virtual inertia of changing the radius of 


gyration of the model from 1.26 feet to 1.17 feet. 


4.2.2 Damping 

The damping was computed using the same analysis program as for 
heave with the appropriate constants changed, It was nondimensionalized 
by multiplying by the square root of the length times the value of grav-~ 
ity and dividing this square root by the displacement times the length 
squared, The values of damping bs a function of Froude number and fre=- 
quency are shown in Table VIII. The damping is plotted in Figure XI to 
show its variation with Froude number, Figure XII is a plot of the zero 


speed conditions, also showing the effect of the reduced radius cf gyra= 


tLON. 


® 


4.3 Theoretical Results 
Based on the coefficients of transformation as presented in Table IV, 

Porter's theory was applied to determine the coefficients of damping and 
virtual mass in both pitch and heave, These coefficients were determined 
at nine frequencies for each station, and the results were numerically 
integrated by strip theory to obtain the coefficients of virtual mass and 
damping for the model. These results are presented in Table IX. These 
Same results are plotted in Figures XIII-XVI to compare their values with 


the experimentally obtained zero speed values, 
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(1) Froude Nunber 





W 
4.13) 
4.909 
6.02 
7210 
7,85) 
Beiter 
9,817 

11.220 
12,083 
13.090 


(2) Froude Number 
W 
4 13h 
4.909 
6.02 
710 
7.85) 
8.727 
FOL 
11.220 
12.083 


13.090 





063 
A 


ol AC Ene 
1.712 
1,82) 
1.788 
eee 
1.718 
1.702 
1.759 
1.80 
1.891 


Table Vv 


a eed 


ove 


1,86 
9250 
8.65 
6.148 
Bee 
4.27 
4.18 
3226 
7208 
8,11 





Max. % Error 





(3) Froude Number 
ih 
413k 
h.909 
6.02 
7e1h0 
7.85) 
6.727 
9,817 
11.220 
12.083 
13.090 


(4) Froude Nunber 
Ld 
4134 
4.909 
6,02 
7Lh0 
7.85) 
8.727 
9,817 
ace 0 
12,083 


13.090 


Table V (continued) 


Oot 
A 
3.76 
2.692 
2,282 
2,113 
1.978 
1.845 
et 
uP ls 
1.75k 
1,862 


0.0 (Parallel) 


A 
3.876 
2.85 
2ec9 
2.186 
2.05) 
1.923 
1.870 
1.905 
1.879 
Less 


a2 B= 





Table V (continued) 


(5) Froude Number 0.0 (Perpendicular) 


w A Al Max. % Error 


413k 2.923 2.220 6.88 
4.909 2.696 . 2.02 7.03 
6,02 2.267 elo 4.10 
7.140 2.097 ; 1.590 3.39 
7.85) 2.030 1.539 4,28 
8.727 1.958 - 1.48) 4.75 
9.817 1.900 | 1.440 5.21 
11,220 1.901 1 bbl 6.8), 
12.083 1.9149 1.478 6.57 
13,090 2.077 Abie 3.90 
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Table vi 





=305 


DAMPING COEFFICIENT IN HEAVE 


(1) Froude Number 3 
v B BS Max. % Error 
4.13) Ho Sie 2.253 Bech 
4.909 6,180 1.900 1.38 
6.02 5.90 : 1.841 4.62 
7 Lh0 6.49 1.99h 2.37 
7.85h 6.637 2.039 3.92 
8.727 6.655 2.0L5 3.61 
9,817 6.120 1.973 Ble 
A220 5.929 1,821 4.00 
12.083 Osrio 1.897 11.12 
13.090 7295 2.2h1 ‘9.27 
(2) Froude Number 2 
w B BY Max, @ Error 
413 4.82 1.488 10.78 
4.909 54279 1.621 15.48 
6.02 4.910 1.509 15.32 
7.140 5.913 1.828 11,88 
785k 6.213 1.910 7.62 
8.727 6.281 1.931 4,00 
9.817 6,182 1.901 Ged 
11.220 5 2636 he f/Sal 46.20 
12,083 5.138 1.578 49.70 
13.090 87h 1.199 100 .00 





(3) Froude Number 





Ww 
413k 
4.909 
6.02 
7210 
7 85h, 
8.727 
9.817 

5220 
12.083 
13.090 


(4) Froude Number 


W 


413k 
4.909 
6.02 
710 
7.85 
8.727 
9.817 
11,220 
12.083 
13.090 





Table VI (continued) 


ai 
B 
5.415 
7.302 
6.129 
6.832 
7.063 
7.05 
6.463 
5.28 
h. 305 
4.493 


QO (Parallel ) 


B 
5.911 
8.119 
6.660 
Teo 
12390 
76034, 
6.112 
5.550 
4.991 
3.431 


one 








(S) Froude Number 


W 
4.13 
4.909 
6,02 
7.10 
7.85) 
O,f27 
9,817 
11,220 
12,083 


13.090 


Table VI (contimed) 


QO (Perpendicular) 


B 
S25 
7430 
8,461 
8.30h 
8,22 
8.005 
(Re 
5.806 
Ne (ee 
36773 


4 


a30= 


Max, % Error 


10.13 
10,52 
7.97 
1213 
Ware 
Mage 
10,86 
22.10 
26.52 
45.30 


VIRTUAL INERTIA COEFFICIENT IN PITCH 


ae SSE ae ee ee ew a SS EE 


(1) Froude Number 
w 
413k 
4.909 
6.02 
7sa0 
7.85) 
8.727 
9.817 
11,220 
12.083 


13,090 


(2) Froude Nunber 
w 
4.13) 
4.909 
6,02 
7 uO 
7.85) 
8.727 
9,817 
my220 
12,083 


13.090 


2.986 
2.528 
2.598 
2,667 
2,605 
2.558 
2.593 
2.h3h 
2.55) 
2.403 


Table VII 


=33= 


Max. 3 Error 


aan 
6.40 
3.15 
03 
T, 36 
1,68 
bis 
Bi (ak 
4,32 
6, 32 


Max. # Error 


vada 
Lehi2 
280 
3G Bs 
Su, ila 

ATi. 

AIC 

ihe, 
4.67 
5.86 





(3) Froude Number 
W 
4.13 
4.909 
6.02 
7 a0 
7.85), 
Op 727 
9,817 
11,220 
12.083 
13.090 


Table VII (continued) 


a 
A 
4,020 
3.538 
3,019 
2.715 
ra (Salk 
2,12 
Zane 
2.200 
3 She 


ZO 


(4) Froude Number O (Parallel) 


W 
4.13) 
4.909 
6.02 
7.140 
7.85) 
Omer 
9,817 

11.220 
12.083 


13.090 


A 
36023 
3.591 
BaLtOT 
2,821 
2,656 
2.519 
2,118 
agen? 
ALY, 
2.19 


aalives 


4.65 
Bpoh 
4s 
ec bil 
2.70 
1.02 
6,09 
7.31 
9.35 
9.81 


Max, % Error 
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(5) Froude Nunber 
Ww 
4.134 
4.909 
6,02 
7.10 
7.85) 
Gaile? 
9,817 
11.220 
12.083 
13.090 


(6) Froude Number 
ka 
4.13 
4.909 
6,02 
7.1h0 
7.85) 
8.727 
9.817 
11,220 
12.083 
13.090 


Table VII (continued) 


QO (Perpendicular) 


A 
399 
3.296 
BROLy 
2.622 
2.1480 
21,08 
2.491 
2.273 
2535 


2ap>3 


Q (Reduced R) 


[> 


3.079 
3.062 
2.67) 
24357 
25 BSS 
2.13 
2,07 
2.037 
1.998 
2,049 


Be 


Max. % Error 


3-10 
1.40 
ec 
ot 
1.59 
Tea fill 
2.69 
4.53 
3.67 
i aGal 


Max. b Error 


029 





Table VIII 


DAMPING COEFFICIENT IN PITCH 


(1) Froude Number ,3 


W B B! 


W B 2 Max. # Error 
4.13) 6.651 073 32.30 
4.909 8.514 09) 114530 
6.02 6.157 .068 13.95 
7.10 5,887 : .065 8.92 
7.85) 5.876 .065 4.95 
8.727 mis ; .057 2ai20 
9.817 5.130 Oy 2.50 

11.220 4.343 048 52.60 
12,083 25020, 2022 100.00 
13.090 4.762 OBE: 40.50 
(2) Froude Number 2 

w B BY Max, @ Error 
413k 8.52) 09), 4.42 
4.909 7.82) 086 2.58 
6.02 5.468 ~060 6.70 
7.10 5.863 2065 Baa 
7.85) 5823 06) 86 
on (all 5.189 2057 1,00 
9.817 4.981 2055 3.148 

11.220 3.588 03 13,60 
12.083 2,921 032 21.60 
13.090 2.003 022 100.00 


wae 





(3) Froude Number 





W 
4.134 
4.909 
6,02 
7.1h0 
7.85) 
Geter 
9,817 

11,220 
12,083 
13.090 


(4) Froude Number 


W 
4.134 
4.909 
6.02 
7.140 
7.85) 
8.727 
9.817 

11.220 
12.083 


13.090 





Table VIIT (continued) 


0.2 
B 


3.087 
TPS Ke) 
5 abl 
6,0),8 
6,318 
5.900 
5 lish 
4665 
3320 
1.276 


O (Parallel) 


B 

820 
3.99k 
5.508 
6.173 
6,091 
5.637 
5.102 
4.205 
aya 
1.491 


=~ 3/= 


Max, % Error 


14.95 
1.01 
1.86 
2.48 
2 
3.72 
5.81 
9,60 

13,36 

20.82 


Max. a Brror 


ns 
aah 
02 
o53 

o 38 
55 
2.66 
Sale 
2 140 
14.40 





Table VIII (continued) 


(5) Froude Number 0 (Perpendicular) 


wW B Bt Max, % Error 
4.13 ~2.723 ~.030 1.40 
4.909 1,351 O15 1h. 36 
6.02 6.197 068 2.57 
fia luale Nees .081 1.0 
7.85), 7.080 .078 10 
eager Ceooe .070 88 
9.817 5.552 - ~060 1.15 
11.220 7.145 079 96h 
i? .083 3.53h «039 10.83 
13.090 10.385 111 23.60 


(6) Froude Number 


w B B' Max. % Error 
4.13 1.87 020 16,83 
4.909 3.886 043 5.90 
6.02 6.085 .067 2.26 
7.140 6.550 2072 1.43 
7.85) 6.155 eral 1.25 
8.727 Only .068 a 
9.817 5.566 061 oll 

11.220 5.395 059 IZ6 
12,083 4.701 .052 5.10 
13.090 4.578 050 15 


O (Reduced R) 


~38- 














Frequency 


4.11 
5.82 
7013 
8,23 
9.20 
10,08 
10.89 
11.6) 
12,80 


Virtual 


Mass 
2,823 
22h 
2.295 
aay 
2.260 
2208 
2.310 
2430 
26369 


Table IX 


THEORETICAL VIRTUAL MASS AND DAMPING 
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FUNCTION OF FROUDE NUMBER 


FIGURE V 


IN HRAVE AS A 
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VIRTUAL MASS IN HEAVE 
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FIGURE VI 


O°S 


O°€T O°CT 


<3 


| (aeTnoTpuedseg Oiesray! || | 1 
(TeTleteg | 


A2uenbeiy Teuorysuswtp—uoy 
O°" Oe 0°? 
£puenbeidy uvewmey 
OE) OFFER 20 1 Oe, O8 BONG OFS O'F 





é 


0 ‘ga 


CET 


O°t 


o 
7 
Ee 
oO 
a-- 
A 
oe 
Oo | 
= 
& 
es 
«{ 
pte - 
| 
= 
@ 
= 
: 
a ' 





ran 
SSey TUN ITA 


° 
ow 


at 
e 
cw 


i 





FUNCTION OF FROUDE NUMBER — 


FIGURE VII 


AS’ A 


DAMPING IN HEAVE 
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FIGURE Vit 
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DAMPING IN HEAVE AT ZERO FROUDE NUMBER 
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FIGURE IX 


VIRTUAL INERTIA IN PITCH AS A FUNCTION OF FROUDE NUMBER 
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FIGURE X 


VIRTUAL INERTIA IN PITCH AT ZERO FROUDE NUMBER 
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FIGURE XI 


FUNCTION OF: FROUDE NUMBER 


AS A 


DAMPING IN’ PITCH 
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FIGURE XII 


DAMPING IN PITCH AT ZERO FROUDE NUMBER 
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FIGURE XIII 


PEED VIRTUAL MASS IN HRAVE | 


COMPARISON OF THEORETICAL AND EXPERIMENTAL 
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FIGURE XIV. 


NG IN HRAVE 


PEED DAMPIN 
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FIGURE XV 


RTIA IN PITCH — 
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COMPARISON OF THEORETICAL AND EXPERIMENTAL 
ZERO SPEED 
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FIGURE XVI 
ZERO SPEED DAMPING IN PITCH 


COMPARISON OF THEORETICAL AND EXPERIMENTAL 
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Chapter V 


DISCUSSION OF RESULTS 


5.1 Experimental Errors 

The resultant curves of virtual mass and damping as shown in Fig- 
ures V - XII show the respective values for the Mariner, In comparing 
these curves with those of other Similar type large hullforms, these 
curves show the same trends in their values, These curves are gener- 
ally believed to give the values for the Mariner at the various speeds 
for the given excitations, However, there are several experimental 
errors or limitations which may limit their accuracy. 

Initially the magnet was supported by a wire. When the tests 
were conducted, the random vibrations were so great as to almost ob~ 
scure the actual record of the model's motion, For this reason stiff 
supporting members were used, This greatly reduced the vibrations, but 
the inherent vibrations in the towing carriage structure and members 
could not be eliminated. These vibrations became greater at the higher 
speeds, They were also affected by the rate at which the model was 
brought up to speed, being much greater when the model was quickly ace 
celerated to the desired speed, In general the vibrations were high 
frequency with little effect on the results, but since the vibrations 
were so random, it is likely that there were enough low frequency come= 
ponents to affect the results, 

Another source of error was introduced by the procedure used in 
recording the results, Initially the response was recorded at 3.75 
feet per second, At this low speed a one-kilocycle-per-second square 
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wave would not reproduce on the tape. In transferring the response to 
a tape at 30 feet per second, it was observed that noise, primarily 
sixty cycles per second, was picked UPe The sixty cycle noise should 
have had no effect on the results, but any lower frequency components 
would have. This effect was considered to be very minor, 

As mentioned earlier, : second circuit was constructed to record 
the actual instant at which the model began to drop, as evidenced by 
the magnet parting from the metal plate, By digitizing this response, 
it was hoped to be able to pinpoint the proper point on the model ree 
sponse data and not ee any errors by starting the analysis at 
the wrong instant. The system was generally satisfactory, but it was 
sometimes difficult to be sure if the chosen starting point was within 
two or three thousandths of a second of the actual one. To find the 
effect of missing the proper starting point, some of the data was ana- 
lyzed at points away from the initial point, The virtual mass coef-= 
ficient showed little sensitivity to small changes, while damping 
tended to be much more dependent on the point at which analysis was be<« 
gun. Thus, this is a possible cause for the errors being generally 
larger in damping than in virtual mass. 

Because of the use of the rotary pitch bearing, friction was nege 
ligible, In heave a definite friction factor was present, However, as 
discussed in Appendix A, the friction had no effect on the response 
curve and thus virtually no effect on the final values of virtual mass 
and damping, To increase the validity of the results, this friction 


factor should have been included in this analysis program. 
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5.2 Character of the Results 

The resultant curves of virtual mass and damping show the general 
form as expected, As predicted by Grim, the curves of virtual mass 
should drop sharply from infinity to a minimum and then rise to an ase 
symptotic value. The infinite point should occur at approximately a 
radian frequency of On k, While the range of frequencies investi- 
gated do not show this point, they do definitely show this trend as it 
is seen that the curves tend toward infinity at frequencies inversely 
proportional to speed, the higher speeds rising at lower frequencies, 

The values of dsp tiie aNeowSiamed the same trend as for other 
ships. Theory predicts that damping should go to zero at higher free 
quencies, However, experimental techniques have shown damping to rise 
at higher frequencies, This indicates that the theory is only good at 
low frequencies, In some cases the damping appeared to still be 
falling to zero in these experiments, However, the higher speeds 
showed a definite rise, The lower speeds were not tested at higher 
frequencies, as this is the general region of towing carriage equip= 
ment natural frequencies. In this region the validity of the results 
was somewhat questionable anyway. 

For the curves in general a study of the results shows the agree- 
ment of the results in the vicinity of the natural frequency of the 
model to be very good, At frequencies removed from this value the er- 
rors went up considerably. This indicated that this method is accurate 
only near the model's natural frequency. Also, as mentioned earlier, 


the errors in virtual mass were much less than for damping, indicating 
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that the procedure may be better for determination of added mass than 
for damping. At the higher frequencies this effect became very pro- 
nounced but was also affected by equipment natural frequencies, 

Kerwin and Narita noted a sudden drop in damping at approximately 
eleven radians per second, This was attributed to the natural frequency 
of the heave rod, In this series of tests a much stiffer heave rod was 
used in an effort to extend the range of frequencies, This sharp dip 
appears to have been eliminated, but there is enough fluctuation at free 
quencies above eleven radians per second to render these results of 
little value, Thus this indieaees an even stiffer heave rod would be 
desirable, 

To overcome the problem of frequency range limitations, there is 
also the possibility of altering the model's natural frequency. This 
was investigated by changing the position of the ballast weights to 
change the radius of gyration. In analyzing the effect of this alter- 
ation, the damping did not appear to change but fell somewhat between 
the two other zero speed conditions. The virtual inertia did change as 
expected and is shown in Figure X. 

The results obtained in pitch did not depend on whether the model 
was initially bow up or stern up. In heave the results were all ob= 
tained with the model lifted up. If a method were used to put pressure 
on the model in order to start with it initially down in the water, it 
would probably produce the same results as when it started lifted up. 
This supports the principle of linearity which was assumed in this ene 


tire procedure. 
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5.3 Theoretical Results 
5.3.1 Transformation Coefficients 

The values of the transformation coefficients ao,4) were found 
which mapped a cylinder into the — of the Mariner, The agreee 
ment between the mapped sections and the actual sections was very good 
in the vicinity of amidships, while the greatest excursions were near 
the ends. For example, at station nine the maximum difference between 
the mapped section and the actual section was only 0.08 feet with an 
average excursion of 0.039 feet. The worst fit occured at station one, 
which had a maximum excursion of 0.9315 feet and an average excursion of 
0.422 feet, The average excursion for the entire ship was 0.168 feet, 
A plot showing the difference in the mapped and actual fom at station 


one is given in Appendix C, 
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From Figures XIII and XIV one can see that the theoretical values of 
virtual mass and damping in heave compare quite favorably with the exe 
perimentally derived values, The virtual mass:values show closer agree= 
ment than do the damping values. However, the difference between the 
theoretical and experimental damping is of ane same order as in the two 
experimental values found from essentially identical response curves, 

This shows the extreme sensitivity of damping to any minor differences, 
and perhaps strip theory is not a precise enough method to evaluate 
damping, unless very fire station spacing were used, 

In pitch the agreement between the experimental and theoretical 


values of virtual inertia was better than the agreement in damping, 
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However, in both cases the variations were greater than for the corres- 
ponding values in heave. This variation was believed to be the result of 
using strip theory on the Mariner hullform. The Mariner has a bulbous 
bow and a wide flared stern, These seeene produced a comparatively 
large value for both added mass and damping. To use strip theory in 
pitch, each section is multiplied by the distance from the center of 
gravity squared and then integrated over the length of the model. These 
extreme bow and stern sections had a very large arm to the center of 
gravity and thus produced very great factors. In the actual case there 
are two reasons why the strip theory is questioned. Near the bow and 
stern the flow is not two-dimensional, as it is near amidship, and thus 
the use of strip theory is not as valid, For example, near the ends the 
water can vent itself not only at the sides, but also at the ends, which 
Should reduce the values found there, These factors may have caused the 
theoretical values to be higher than the values found experimentally. 

As in heave, it is believed that much finer spacing of sections 
should be used, especially near the bow and stern in order to obtain bet-~ 
ter results for the theoretical values. This will not correct the error 
resulting from three dimensional effects, but would reduce errors caused 
by too large an integrating interval in the vicinity of these critical 


regions, 
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Chapter VI 


CONCLUSIONS 


6.1 Evaluation of Step Method 

It is considered from the results of this work that the step dis- 
turbance method is a good method for obtaining added mass and damping 
coefficients, It has several advantages in that the experiment is very 
easy to set up. By use of a device, such as a magnet, it requires only 
low potential devices for power requirements, which are readily avail- 
able. The calculations involved would be virtually impossible if done 
by hand, but they are readily adaptable to computer operation, One 
problem that may arise with some models is the problem of wall reflece 
tions in a narrow tank, This was not a factor in the Mariner tests, 
but by orienting the model to the centerline of the tank at various 
angles, this problem can be overcome, 

It is concluded further that the process is better suited for 
added mass than for damping. This is a factor of the experiment and 
the analysis program, Damping variations were much greater because of 
vibrations or not having the exact starting point in the data. From 
the analysis program, it was seen even with exact generated data, damp= 
ing showed variation from the actual values, However, if care is taken 
to get minimum vibrations and if the data is sampled at fine enough in- 
tervals, it is believed that the process is valid for damping as well 
as for added mass. 

A limiting feature of the experiment is the dependence of the re= 
sults on the model's natural frequency. This, coupled with the 
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interference of the supporting structure's natural frequency, could 
limit the range of reliable results, By using a very stiff heave rod 
and bracing it, its natural frequency could be made high enough to not 


interfere with the results, 


6.2 Experimental Results 

It is observed from the experiments that both added mass and 
damping are speed dependent. However, in the vicinity of the natural 
frequency of the model, this dependency is very slight. The speed dee 
pendence was generally greatest at the lower frequencies, 

Virtual mass was affected = model's natural frequency, while 
damping showed little dependence, Thus it appears possible to extend 
the useful range of frequencies in damping by testing the model at vari- 
ous natural frequencies, For virtual mass calculations the proper 


natural frequency desired must be used for the model, 


6.3 Theoretical Conclusions 

The method of mapping a cylinder into a ship section is very good 
if sufficient coefficients are used, With the exception of highly ir 
regular sections, almost perfect agreement between mapped and actual sec~ 
tions can be obtained, However, the extent of the calculations make it 
almost impossible to perform them by hand, but they are readily adapt 
able to computer operation. 

The virtual mass and damping coefficients derived from this method 
were believed to be accurate, as evidenced by the results obtained in 


heave. However, to apply strip theory to the results in pitch, it is 
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recommended that finer spacing be used at the bow and stern, especially 
for irregular sections such as are on the Mariner, By this method a 
more accurate representation of the virtual inertia and damping coef- 
ficients at the ends of the model could be obtained, Thus, when multi- 
plied by the distance to the center of gravity squared, these sections 
would probably not produce the effect of having the bow and stern domi-e 
nate in obtaining the coefficients of virtual inertia and damping for 


the model, ; 
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Chapter VIT 


RECOMMENDATIONS 


7el Experimental Technique 

The major problem faced in analyzing the data was obtaining the 
data in digital form, The analog to digital program used was very ef~- 
ficient, but there were Several mistakes made by the operators which, 
if not made, would have greatly aided the obtaining of the digital 
output. “ 

The first problem was recording data at too slow a speed and not 
being able to record the digitizing signal, The tape recorder has to 
run at a speed of thirty feet per second to record a one~kilocycle-per- 
second Square wave properly. 

Each run took approximately twenty-two seconds to digitize, If 
data were recorded with this spacing between runs, the digitizing cone 
sole should be set up to run continuously, and there would be no lost 
time in setting the start'of each run when the computer was ready for 
more data, 

In setting up the digitizing console, one has to set up an ame 
plifying and an offset section so as to keep the data between plus and 
minus eight volts. It is best to get as wide a spread as possible in 
order to get the maximum resolution of the data, In order to get this 
spread and not have any data be clipped, it is essential that all runs 
be recorded to the same scale, Thus in setting up the recording equip= 
ment, the zero position and maximum deflection position should be made 
the same on the recording equipment for every run, Otherwise reduced 
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scales must be used to keep all data within the allowable range, with 
some loss in accuracy, 

The final point to be made about digitizing concerns the treat- 
ment of the magnet parting signal. “the method employed by the opera- 
tors in this instance produced several problems that could have been 
avoided. The first problem was that in putting the one kilocycle 
Square wave on the tape, the instant of magnet parting was missed in 
some cases as the Square wave was put on the tape too late. These 
runs had to be digitized a second time to get useful data. This was a 
waste of computer time that was necessary to get the maximum amount of 
data analyzed. Finally, in scanning the record of the parting signal, 
it is possible that the starting point may have been missed by a couple 
thousandths of a second, ‘this Same accuracy could have been obtained 
by simpler methods, The best method appears to be to wire the square 
wave generator between the recorder and the magnet such that when the 
magnet parts, the square wave begins recording. This process should 
have negligible time delay. 

The other big problem was that of random vibrations in the data, 
As stiff a supporting member as possible should be used for supporting 
the ship and supporting the magnet. This would not only reduce the 
vibrations, but also would not limit results because of the equipment's 
natural frequency. Also, the model should be accelerated up to speed 
as slowly as possible to minimize the vibrations caused by this accel- 
eration, 


To increase the validity of the results found in heave, a device 
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employing rotary joints and stiff members may have the benefit of less 
friction than the heave rod, The inertia of the mechanism may be dif- 
ficult to determine, but it may be a better mechanism for this type of 


experiment, 


72 Future Studies 

The method employed here is considered a good one and could be 
applied to further studies in ship motions, It could possibly be used 
to obtain the values of the coupling coefficients between pitch and 
heave, This could be done with a similar arrangement in which the model 
were allowed both types of freedom, and the two motions recorded, The 
method could also be extended to further studies of damping. ae could 
be done by measuring the relation between wave amplitude at the sae 
and the amplitude of the model's oscillation in heave, The same linear- 
ity assumptions could be employed, and the results compared with those 


obtained by forced oscillation tests. 
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APPENDIX A 
DISCUSSION OF EQUIPMENT FRICTION 


In choosing the components of the towing apparatus that allowed 
pitch and heave motions of the ship, and in choosing the measuring dee 
vices, it was desired to eliminate or reduce friction to a minimum, 
Before any connection or measuring pians were made, this was the first 
consideration, | 

The problem in pitch was not very great. The only moving cone 
nection, which also contained the measuring device, was the pitch 
bearing, It is a small bearing rotary device and contributed no noe 
ticeable friction. 

In heave the problem was much greater, The first idea was a 
mechanism employing straight bars connected by small bearing pivots, 
While there would be several such pivots, each would have essentially no 
friction, and the net effect would approach frictionless operation, The 
biggest problem was the system complexity required to insure heave motion 
only. Also, there was a problem of determining the exact system inertia 
in motion, 

Because of the above problems, it was decided to investigate the 
use of a heave rod, Previous heave rods employed had a large amount of 
friction, A new one was constructed using a square rod moving on a 
system of small roller bearings, Friction tests were conducted on it by 
a weight-dropping method, Various weights were attached to the rod, and 


the acceleration was measured as the rod fell. Comparing this acceleration 
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to one gravity unit, one could determine the friction in the mechanism, 
Furthermore, the weights were hung at’ fifteen and thirty degrees to the 
vertical to give the effect of the horizontal force induced by the model 
going through the water at various speeds, The results of these tests 


are tabulated in Table X, 


Table X 


HEAVE ROD FRICTION 


Angle of Weight . _ Average Friction 
Oh »385 lbs, 
15° .506 lbs, 
30° 822 lbs, 


To gain insight as to the effect of this friction, a test was con~ 
ducted at zero speed with the heave rod removed, The model was attached 
to the towing carriage such that it was free to oscilldte on the linear-= 
syn only, It was set to the same amplitude as was used in the heave 
tests, When the model was released, and the record of the heave com- 
pared to that with the heave rod, there was no apparent difference in 
the response, Accordingly, the effect of heave rod friction was nege 
lected in the final calculations, 

As this test was conducted at zero speed and as the weight drop 
tests showed considerable increase of friction with a horizontal com= 
ponent, this may have been a source of slight error at higher Froude 
mimbers, However, this effect was considered secondary and was nege 


lected, 
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APPENDIX B 
EFFECTS OF WALL REFLECTIONS 


A question raised early in the testing was how valid the results 
would be or whether the effects of wave reflections in a narrow tank 
would be such as to render the response meaningless, Because of this 
possibility, two procedures were available to get better results at 
zero speed. One was to test the model at different orientations in the 
tank, and if this proved ee to make the tests in a larger 
medium, such as a swimming pool, 

The model was first tested parallel to the tank. The heave ree 
sponse obtained is shown in Figure XVII. The model was then tested at 
an orientation perpendicular to the tank, and the response is shown in 
Figure XVIII, Similar results were obtained for pitch, In both cases it 
is noted that the model oscillations damped to zero before any wall ree 
flections affected the motion of the model, It is true that the exe 
tremely long waves returned to the model almost instantly, but their 
effect is negligible, and it is the wave of the same order of magnitude 
as the ship length with which we are primarily concerned. 

It is interesting to note that there is a considerably longer 
zero region when the model was placed across the tank than when it was 
oriented parallel to the tank. This indicates that wall reflections 
must be considered in this type of analysis, but fortunately, for the 
Mariner model the oscillations completely damped out prior to wall ree 
flections, Thus the results of the analysis were essentially the same 
for each condition. 
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APPENDIX C 


DETERMINATION OF COEFFICIENTS OF TRANSFORMATION 


The determination of the coefficients of transformation of the 
semicircle in one complex plane to the ship section in the second com- 
plex plane proved to be very adaptable to computer operation, The IBM 
709, computer at the Massachusetts Institute of Technology was used 
for this purpose in this project. In general, the mapped cylinder gave 
a very close approximation to the actual ship section. The method used 
consisted of selecting tentative values of a}, 83, and ac by reference 
to Landweber and Macagno (9). The computer program was then entered 
with these values and offsets of the actual ship section, An iteration 
procedure written into the program generated new values of the coef- 
ficients and offsets which were compared with the actual ship section. 
The iterations were continued until a reasonably small least mean 
square error between calculated and actual offsets was obtained, The 
maximum deviation used in this project occured at station one, Fig- 
ure XX in this appendix compares the actual station one with station 


one derived from the transformation coefficients, 
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FIGURE XIX 


COMPARISON OF THE ACTUAL MARINE? STATION ONE 
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After the main body of work on this project was completed, a new 
method for determining the transformation coefficients, 2on¢]>» Was de- 
vised, This method, together with the new values of the transformation 
coefficients for the statzeneror the Mariner hullform, are presented in 
this appendix, 

Equation (5) represents the conformal transformation of the unit 


circle to the section in question. In equation (5) 


« 


C ea dys ae . (17) 


Z=xX +iy = iretoe Gis) 


In this work the y axis points down, and the angles 9 and 8 are measured 
from the vertical, X is the offset of the actual form, and y is the 
waterline, Therefore, when 6 = 8 = 90°, x is the half-beam of the sec# 
tion, and y = 0, When 6=6 = O°, x = 0, and y is the draft of the sec~ 
eon. 

By substituting equations (18) and (17) into equation (5) and sepa- 


rating the real and imaginary parts, the following equations are obtained: 


N 
Sane ain n “2n+1 sin(2n +1 1 
X= Poll eta, (-1) ase ( )8 ( 9) 


n=O 
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a 
\ ( =) ala ene: 


y= pcose +? ens cos(2n + 1)6e (20) 


It is then necessary to solve for @ and Aon+1 Which satisfy these equa- 
tions, 

If M points on a section outline are available, there are M pairs 
of equations of the above form, or a total of 2M equations. The values 
Of ao,41 must be the same in ail of these equations, but the values of 

will be unique for each given point. 

If the half-beam and draft of the section are known, an additional 


equation is available, As mentioned above, at © = 90°, 





a 
\ en+1 
half-beam = b =p + (21 ) 
/ eOnt+L 
n= P i 
At © = 0°, 
N a 
draft =H =p + (Gi a a (22) 
n=0 P 


These results lead to the equation 


| 


eee b 
pl -F) + (1 +2) +30 - P+ ee. 0 (23) 


We now have 2M+] equations which mst be satisfied. 
In the new method of solving for the coefficients a5,,,, the co- 


ordinates of the M points, as well as the half~beam and draft of the 
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section, are inputs to the computer program, Three alternatives are now 
available, If estimated values of the coefficients ao,,4) are available, 
they may be entered. In lieu of this, if estimates of the angles 6 are 
available, they may be entered as inputs to the program. If neither the 
coefficients nor the angles 6 are available, the program will generate 
an estimate of the angles 6. | 

Assuming that either the angles § are inputs or that values of 8 
are internally generated, the program solves the eM+] equations for Pp 
and aoy,4 1 to give a least mean square error solution between the gener- 
ated section and the actual ship section, The values of the coeffic- 
ients are then used with input values of x and y to determine a new 
value of 6 for each point, The program then re-enters the least mean 
Square process and iterates until the desired number of iterations is 
completed, If values of Qon41 are inputs, the program computes values 
of & before entering the mean square process, 

It should be noted that the program requires the number of co- 
efficients desired and the number of iterations desired as inputs, 

As stated above, the program finds values of p. However, since 
we are interested in the case of the unit circle, the program is nor- 
malized with P equal to one. 

Using this new method, new values of aj,,, were determined for the 
Mariner hullform. These coefficients are presented in Tables XI and XII. 

It is believed that this improved approximation to the hullform 
will not appreciably affect the values of damping and virtual mass for 


the Mariner, However, in future work it is recommended that the new 
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method presented above be used, 


A listing of the analysis program in MAD computer language is at- 


tached to this appendix, 


Station 


OO NM & Ww NY fF ©€ 


sgl. 

-.858571 
~ 675965 
~ 479528 
-.286588 
~-.117882 
~.000)35 
+ ,076025 
+,121668 
+1399), 
+.150007 
+,150031 
+.1511h) 
#150825 
+ .1)8566 
+ 139671 
+,1174)10 
+ 069663 
~ 026919 
~ 216518 
-.562587 


+ 227269 


Table XI 


a3 
~ 118183 
-.085363 
~ 04,3086 


= ,022750- 


-,0208L) 
-.035251 
-.05 3853 
-.078578 
~.102216 
- 11900 
- 12,76 
~,122959 
~.110575 
-,085017 
~ 08057 
~.002h,30 
+ ,05 3825 
+ .106876 
ape rloe ya 
+ ,131206 


+ 251236 


A-1l 


TRANSFORMATION COBFFICTENTS 


a5 
~,018627 
~ 51353) 
+ ,01033 
+ ,019):79 
+ ,019821 
+ 017899 
+ ,016527 
4.01077 
+ 001109 
-,001388 
- 004275 
~ 003435 
~ 001117 
+ ,0055 37 
+ 016075 
+ 02435 
+ 032166 
+ 01,3512 
+ 074115 
+ ,0982))5 


=— 12 O12 


a7 
+ ,006032 
+ ,006935 
+ ,003225 
+ ,003222 
+ ,005278 
+ 002628 
+,001h57 
-,00067) 
-,001911 
~ ,002853 
-.00)280 
-,00L10) 
- ,00)822 
-.00,659 
~ ,00566 
-,001183 
+, 00002 
~.001900 
+ ,012576 
+ 047890 


+ 048271 


2 
-,000315 
+ ,005787 
+,013838 
+ .016652 
+ ,010582 
+ ,0078h 
+ 005735 
4.002915 


+ ,001 686 


+ 001582 
+, 006575 
+ ,009 33h, 
+.006618 
= ,003023 
+ ,016198 


-.109952 





Table XII 


COMPARISON OF ERRORS WITH THOSE FOUND BY PROCESS 


ia eee oe eee 


a 


A=12 


Average Maximum Average Maximum 
Station Excursion Excursion Excursion Excursion 
0 0.07 0.11 
1 0,09 Oly 0.32 0.93 
2 0.16 0.33 0.17 0. 32 
3 0,13 0.25 0.20 0.54 
hy 0.12 0.23 One: 0.0 
5 0.10 0,22 0.13 Opec 
6 0,08 0.17 0.08 Oj aly 
il 0,06 0.13 0.06 0.18 
8 0,03 0.07 0.07 0.13 
9 0.02 0.07 0.02 0.08 
10 0.05 0.12 0.08 0.23 
1 0.06 Oza 0.06 O.11 
are 0,03 0.09 0.03 0.09 
13 0,03 0.05 0.0) 0,08 
Lh 0.03 0.09 0.0 0.10 
us 0,09 0.25 0.19 0.40 
16 0,08 0.19 0.2) 0.55 
ae Oel2 Oeer 0.19 0.52 
18 0.06 0,11 0,22 0.46 
19 0.15 0.30 0.29 0.69 
20 Ore 0,13 0.07 0.20 





RCOEFFICIENT DETERMINATION FOR HULL FORMS 

DIMENSION B(625sBDIM) sTHETA(6259TDIM) 9X( 30) 9Y(30) sAGAR( 25) » 
} 1 APRI (15) sBETA(30) sCRIT(20) 9A(15) »WYE (30) 9c X(30) sERRX(30)¢ 
2 ERRY(30) »NAME(12) sDEGREE(6259DDIM) sALPHA( 30) sDUM(22>59U0I1MD) 
INTEGER I sJ9K o9MsNoQeCOUNT sPWReSTGNAL 


VECTOR VALUES BDIM = 291925 
VECTOR VALUES TDIM = 291925 
VECTOR VALUES DDIM = 291925 
VECTOR VALUES DIMD = 2591 91D 


START READ FORMAT TITL1sNAME(1)eeeNAME (12) 
VECTOR VALUES TITL1 = $(12C6)*5% 
PRINT FORMAT TITL29NAME(1)¢eeNAME(12) 
VECTOR VALUES TITL2 = $(1H1912C6)*4% 
READ FORMAT INPUT» HoeMeNsQeBEAM,SIGNAL 
VECTOR VALUES INPUT = S(F102693129F10¢69I12)*S 


BDIM(2) = N+2 
DIMD(2) = N+l 
TDIM(2) = Q+1 
DDIM(2) = Q+]1 


READ FORMAT OFFSET» X(1l)eeeX(M) 
VECTOR VALUES OFFSET = $(5F1004)*$ 
READ FORMAT Whose Yll)eeeY(M) 
VECTOR VALUES WL = $(5F1024)*$ 
KAY = BEAM/H 
El = 1e0E-96 
WHENEVER SIGNAL eEe O 
READ FORMAT ANGLE» EX(1)-eeeEX(M) 
THROUGH BEEl»s FOR-I = lols IeGeM 
Bee l DEGREE(Is1) = EX(T) 
VECTOR VALUES ANGLE = $(5F10e¢5) *$ 
THROUGH TWENTYs FOR I = lololeGeM 
TWENTY THETA(I 91) = DEGREE( 1 91)/572295780 
TRANSFER TO TWELVE 
OR WHENEVER SIGNAL eEe 1 
READ FORMAT AV/fLs All)eeeAI(N) 
VECTOR VALUES AVAL = $(5494E14%e6)*& 
THROUGH THIRTYs FOR I = Lle2evleGel2*N—-1 ) 
THIRTY A(2*I-1) = ACT) 
DEN = le + All) 
THROUGH FORTYe¢ FOR I = 32 aT eGel2*Nn-l ) 
FORTY DEN = DEN + A(T) 
RO = BEAM/DEN 
THROUGH FIFTY» FOR I. 
FIFTY APRI( I) = ACT) #RO 
TRANSFER TO TWELVE 
OTHERWISE 
CONTINUE 
END OF CONDITIONAL 
R 
R 
THROUGH INITL»® FOR I = LlololeGeM 
THETA(CI 91) = ATANe(X(IT)/YCI)) 
INITL WHENEVER Y(I) elEe Els THETACI91) = 1457080 
TWELVE THROUGH ONEs FOR J = loleJeGeQ 
WHENEVER SIGNAL eEe 1 eANDe J eEe leo TRANSFER TO Bl 


% 


le2esTleGel2*Nn-l ) 





R 
R 
THROUGH TWO*x FOR I = lols LeGemM 
Gemeety = SINe ( THETA( I oJ)) 
BlIs2) = BlIel) 
B(Is(N+2)) = X(T) 
B(CI+M)91) = COSe(THETA( I 9J)) 
B((Il+M)092) = ~—F((I14+M) 51) 
BistileeMystN+2)) = Y(I) 
THROUGH TWOs FOR K = 39ls KeGel(N+l) 
Bi IeK) = (lle) ePeKI*SING (((2*K)—-3)*THETA( I 9J)) 
TWO BCC I+M) 0K) = —€Cl—-lLedePeK) *#COSe(((2*K)—-3)* THETA I 9J)) 
THROUGH THREEs FOR K = 1lo9l9KeGe(N+1) 
THREE B((2*M+1) 9K) = le + ((—LelePeK)*#KAY 
B( (2*M+1)9(N+2)) = Oe 


EXECUTE GLSQHe(B »sABAR sDUMs (N41) 59(2*M+1) ) 
R 

PRINT RESULTS ABAR(1)eeeABAR(N+1 ) 

R 

RO = ABAR(1) 

THROUGH FOUR»s FOR I = 291s IeGe (N+1) 

FOUR APRI((2#1)-3) = ABAR(I) 

R 
R 

THROUGH FIVE» FOR I = 1lo29 [eGe((2*N)—1) 

FIVE A(T) = APRI(I)/RO 

R 
R 

THROUGH SIXs FOR I = lols IeGeM 

EX(1) = (RO + APRI(1))#*SINe( THETA( I 9J)) 

PWR = O 

DEGREE(IsJ) = 57e¢2958*THETA(I 9J) 

WYE(I) = (RO-APRI(1))#*COSe(THETA(I 9J)) 

THROUGH SIXPRI» FOR K = 3929 KeGell2*N)-1) 

PWR = PWR + 1 


WYE(I) = WYE( I) — (l—-lLledePePWR)#APRI(K)*¥COSS(K*THCTA(C I sJ) ) 
SI XPRI EX(1) = EX( 1) + (lle) ePePWR)*APRI(K)*SING(K*THETA(I 9J) ) 

mre yl) = SAP Se (WYECE) = YOR) ) 
SIX ERRX(I) = eABSe (EXIT )=X(T)) 

R 

R 

PE aA Oj = 0 

BETA(O) = Ov 


THROUGH EIGHT» FOR I = lols IeG,.M 
ALPYA(I) = ALPHA(CI-1) + ERRX(I)#ERRX(1) 
EIGHT BETA(I) = BETA(I-1) + ERRY(CI)*ERRY(TI) 

R 
R 

CRIT(J) = SOQRT.(ALPHA(M) + BETA(M)) 

FXECUTE OUTPT, (J»sROeAsCRIT o9X v9 EX sERRXo YoWYE sERRY sDEGREE o MoN) 

WHENEVER JeEeQs TRANSFER TO FINIS 
R 
R 





MHROUGH CEEl» FOR I = lo2s IeGe(2*N-1) 
-FF] PRINT FORMAT ARCs A(T) 
: VECTOR VALUES ABC = S(1H »£14e6)*$ 
PRINT FORMAT LMS» CRIT(J) 
VECTOR VALUES LMS = $(21H LEAST MEAN SQUARE = 9F14e6)*% 


PRINT COMMENT $0 MEASURED XACT) CALCULATED X(T) X(1) 
ieee ROR MEASURED YC(1) GaSCULEATED  y« 1} bl i 0 DR ne 
20R THETA(I)$ 

THROUGH TENs FOR I = lols IeGeM 

TEN PRINT FORMAT OFVAL » Mite X11) ,ERRX(1T) sYCI)V SWYECI) sERRYC(I) 9 


l DEGREE) 
VECTOR VALUBS@OEVAL = S(1H oF 146955 9F 14069559 F 14069599 
i P14. 6995 9F 1406955 9F 1406955 9F 1406) ¥S 
FUNCTION RETURN 
END OF FUNCTION 


MAD 
EXTERNAL FUNCTION (AsX9BsNoM) 
ENTRY TO GLSQHe 
INTEGER ToJoKoNoM 
THROUGH QQl»> FOR I = lols IeGeN 
X(I) = Oc 
THROUGH QQ1» FOR J = lols JeGeM 
QQl X(T) = X(T) + AlSsT)*ACS9 (N41) ) 
THROUGH QQ2> FOR I = lols IeGeN 
THROUGH QQ2+ FOR J =* lols JeGeN 
BileJ) = Oc 
THROUGH QQ2s FOR K = lslo KeGeM 
QQ2 B(IeJ) = BlIsJ) + AlKelI) *A(K9J) 
THROUGH QQ39 FOR I = lols IeGeN 
WHENEVER B(I9I) eEe Oo 
THROUGH QQ4o FOR K = (I+1)919 KeGeN 
THROUGH QQ5e FCR J = Tels JeGeN 
TEST = BlIoJ) 
BitleJ) = B(K>sJ) 
905 Bi(KeJ) = TEST 
TES? = XCEL) 
Cl he = etn) 
KURee= TEST 
QQ4 WHENEVER B(I51) eNEoe Oes TRANSFER TO QQ6 
END OF CONDITIONAL 
2Q6 STORE = Bl(IoeI) 
THROUGH QQ7s FOR J = Ils JeGeN 
QQ7 Bl ls =meels J) /STORE 
X OD earl So ORE 
THROUGH QG3,5 FOR K = lols KeGeN 
WHENEVER KeEeTeORe Bl Kol) eEe0ee TRANSFER TO QQ3 
STORE = Bi(KolI) 
THROUGH QQ8» FOR J = Ils JeGeN 
B(KeJ) = BlIlsJV*¥STORE = B(KsJ) 
QQ8 WHENEVER Kelels B(KsJ) = =B(KoJ) 
X(K) = X(T)*STORE = X(K) 
WHENEVER KelLels X(K) = —-X(K) 
QQ3 CONTINUE 


FUNCTION RETURN 
END OF FUNCTION 














